Astronomy Sz Astrophysics manuscript no. 091025Coolingv4 


© ESO 2009 


November 26, 2009 





On 
O 
O 

> 
O 

(N 
< 

o 

Oh 
6 

(N 
> 

o 

in 

o\ 
o 

On 
O 



13 



A new radiative cooling curve based on an up to date plasma 

emission code 

K.M. Schure 1 *, D. Kosenko 1 , J.S. Kaastra 1 ' 2 , R. Keppens 1 ' 3 ' 4 , and J. Vink 1 

1 Astronomical Institute, University of Utrecht, Postbus 80000, NL-3508 TA Utrecht; K.M.Schure@phys.uu.nl 

2 SRON, Utrecht, The Netherlands 

3 Centre for Plasma Astrophysics, K.U. Leuven, Belgium 

4 FOM-Institute for Plasma Physics "Rijnhuizen" , Nieuwegein, The Netherlands 



Received . . . / Accepted 



ABSTRACT 



This work presents a new plasma cooling curve that is calculated using the SPEX package. We compare our cooling 
rates to those in previous works, and implement the new cooling function in the grid-adaptive framework 'AMRVAC. 
Contributions to the cooling rate by the individual elements are given, to allow for the creation of cooling curves 
tailored to specific abundance requirements. In some situations, it is important to be able to include radiative losses in 
the hydrodynamics. The enhanced compression ratio can trigger instabilities (such as the Vishniac thin-shell instability) 
that would otherwise be absent. For gas with temperatures below 10 4 K, the cooling time becomes very long and does not 
affect the gas on the timescales that are generally of interest for hydrodynamical simulations of circumstellar plasmas. 
However, above this temperature, a significant fraction of the elements is ionised, and the cooling rate increases by a 
factor 1000 relative to lower temperature plasmas. 
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1. Introduction 

Optically thin plasmas can cool very efficiently, especially 
on timescales relevant to circumstellar evolution. A correct 
description of cooling is therefore important for simulations 
of e.g., stellar winds. The dominant cooling mode for plas- 
mas within the temperature range of 10 4 — 10 7 K is by 
metal line transitions. In many cases, the plasma can be 
assumed to be in collisional ionisation equilibrium (CIE). 
A predefined cooling curve can then be used to account 
for radiative losses, where cooling rates depend only on 
the temperature of the plasma. In the calculation of the 
cooling curve, we can adjust the abundances and try to in- 
corporate as many line transitions as are reliably known, to 
obtain a straightforward and efficient way of treating radia- 
tive cooling. Well-documented cooling curves used in HD 
an d MHD computations are those for example develope d 
bvlDalgarno k McCrayl (11972l)lRavmond k Smith! (|1977l) . 
iMacDonald k Bailevl (ll98lNSutherland k Dopital Jl993), 
and Smith et al.1 (120081). Other approaches, such as that 
of iMellema k LundqvistJ ([2002) , include time- and ion de- 
pendent cooling, with CSM abundances that depend on the 
evolutionary stage of the star, but treat the hydrodynamics 
only in one dimension. 

We use the SPEX package (jKaastra et al.lll996[) to cal- 
culate a new cooling curve, using up-to-date transition lines 
to calculate radiative losses for a plasma. This code is used 
as a spectral analysis code tailored to EUV and X-ray ob- 
servations, in which energy band it is one of the most com- 
plete packages currently available. Because it includes a 
very complete prescription of line emission, cooling rates 
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predicted by SPEX are higher than those of cooling curves 
available until now. 

We implement the new cooling curve, as deduced 
from the SPEX CIE co de, in the AMRVAC frame- 
work (|Keppens et al.| [2003), of which we use the hydrody- 
namics to demonstrate the influence of radiative losses on 
the gas dynamics. As an example, we use this code with 
the new cooling prescription to simulate the circumstel- 
lar medium of a massive star. Plasmas typical of this cir- 
cumstellar medium are at fairly high temperatures and are 
mostly optically thin. Therefore, they can cool very effi- 
ciently. 

The cooling at a certain point in the simulated plasma 
can be calculated from the local temperature and density 
in combination with the prescribed cooling curve. Different 
cooling curves can be chosen, which depend on metallicity 
and/or ionisation fraction (for low temperatures). We pro- 
vide cooling rates for the individual elements, so that the 
cooling curve can be tailored to fit the required abundance 
ratios. 

In Sect. [21 the SPEX code, the abundances, and the 
new cooling curve will be presented and compared to exist- 
ing cooling curves. We show that the SPEX curve causes 
more cooling for temperatures from 10 to 10 6 K, a tem- 
perature range that is very important in the evolution of 
the circumstellar as well as the interstellar medium. 

In Sect. 02 we discuss the range of validity of the CIE 
assumption and discuss the implications of photoionisation 
and shock heating on the cooling rates of the plasma. 

Cooling at lower temperatures cannot be accurately 
calculated with SPEX, because of incomplete data for 
A > 2000 A. In Sect. @J we briefly summarise our 
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treat ment of lower tem p eratur e plasmas, which fol- 
lows iDalgarno fc McCravl (|1972f) . The cooling rates of 
ISmith et al.l (|2008l ), provide a detailed representation of es- 
pecially the low-temperature (T < 10 4 K) cooling of plas- 
mas. We select our low-temperature ionisation fraction to 
ensure the optimal correspondence with their cooling rates. 

In Sect. we show the effect that cooling has on 
the dynamics of the plasma, by using the SPEX cooling 
curve as implemented in the AMRVAC code. Because the 
cooling adds to the resolution requirements of the hydro- 
simulations, the grid-adaptive method used by AMRVAC 
maximises the performance while reaching sufficient reso- 
lution to enhance the higher density contrast and allow for 
instabilities to develop in the cooling shell. In Sect. O we 
summarise the important points and outline the major dif- 
ferences from earlier cooling curves. 



2. SPEX cooling curve 

To generate the cooling curve we u se the SPEX0 package, 
version 2.00.11 ([Kaastra fc Mewdl2000l) . This code is tai- 
lored to fit spectra in the EUV and X-ray regime and can be 
used to calculate the luminosity of a given emission measure 
for plasmas at different temperatures and different choices 
of electron and ion temperatures. For our purpose, we use 
SPEX to calculate the emissivity of a plasma in collisional 
ionisation equilibrium (CIE ). The code is b ased on the 
mekal line-emission model ([Mewe et al.l Il995f) , but many 
updates have since been made. The 15 elements presently 
included are H, He, C, N, O, Ne, Na, Mg, Al, Si, S, Ar, Ca, 
Fe, and Ni. 

We calculate spectra for a grid of temperatures rang- 
ing from 10 3 - 8 to 10 816 K, with steps of logT = 0.04 K. 
Solar abundanc es are assumed, usi n g the abundance ta- 
ble presented in lAnders fc Grevessd ([1989D . which is sum- 
marised here in Table [T] for elements that are taken into 
account. The emissivity Ajy (erg s -1 cm 3 , normalised to 
n e nnV = 1 cm~ 3 ) is calculated by integrating the spec- 
trum over an energy range from 10 _1 eV to 1 MeV. The 
SPEX package is currently one of the most complete spec- 
tral packages, and it is especially complete for temperatures 
higher than - 10 4 86 K (see Sect.©. 

The cooling curve is tabulated in Table [2l and in 
Fig. Q] it is plott ed against the CIE cooling curve that 
was calculated bv lSutherland fc Dopital (|l993l ). Compared 
to their curve, and other CIE cooling curves currently 
available, we find that the cooling efficiency calculated by 
SPEX is higher. We attribute this difference to SPEX tak- 
ing into account more line transitions in both X-ray and 
EUV than p revious studies (5466 l in es in SPEX compared 
to 1683 in ISutherland fc Dopital dl993h ). Additionally, 

1993lT ~ 



though ISutherland fc Dopital (|1993l ) use the same val- 
ues of Solar abundance ratios, they use a slightly dif- 
ferent ionisation balance, which they calculate with their 
MA PPINGS II code. SPEX use s the ionisation balance 
from lArnaud fc Rothenflugl (119851) for mos t eleme nts, apart 
from the one from lArnaud fc Raymond! ([l992) for iron 
(|Mewe et al.lll995l) . A more complete list of Fe L complex 
lines is al so included, based on the HULLAC atomic physics 
package (jKlapisch et al.|[T977h . causing higher cooling effi- 
ciency around 1 0_ K. The differen ces in the treatment of 
Fe are shown in iMewe et al.l (|1995h . 



Additionally, we 
iMacDonald fc Bailevl 



plot the cooling curve by 
(1981) (hereafter also referred 



to as MB81) in Fig. Q] We used this for comparison in 
our hydro-simulations in Sect. El This cooling curve is a 



composite of cooling rates by Raymon d et al 



10 6 < T < 10 s K and by IShapiro fc Moore 



(1976) for 
( 1976 ) for 



10 4 < T < 10 6 K. Below 10 4 K, extrapolation of the cooling 
rates is used, towards zero at 100 K. Their treatment of 
cooling between 10 4 — 10 6 K is a non-equilibrium approach, 
applicable to a low density plasma that is overionised, 
as a result of rapid cooling (t cool < ^combination)- Time- 
dependent cooling, starting from ionisation equilibrium 
at 10 6 K, gives rise to lower cooling rates, as discussed 
in Sect. [31 MB81 base their cooling rates on the inclusion 
of 12 elements (H, He, C, N, O, Ne, Mg, Si, S, Ca, Fe, 
Ni), and use slightly different abundances, which are not 
however responsible for the major differences from our 
cooling rates. 

To apply the cooling curve in the AMRVAC code, we 
renormalise the cooling rates calculated with SPEX by mul- 
tiplying by the electron to proton ratio (ne/nn). We thus 
obtain a value for A n d that we can use in our hydrodynamic 
application: 



A 



ii<i 



-A 



N- 



(1) 



The ratio n Q /nn is a temperature-dependent parameter 
that is also calculated by SPEX and is given in the last 
column of Table [2] The cooling rate of the plasma is given 
by 



L = J nHnnA hd (T)dV (erg s" 1 ). 



(2) 



Table 1. 

(1989) 



Solar abundances from lAnders fc Grevessi 



element 


log(n/n H ) 


element 


log(n/n H ) 


He 


-1.01 


Al 


-5.53 


C 


-3.44 


Si 


-4.45 


N 


-3.95 


S 


-4.79 


O 


-3.07 


Ar 


-5.44 


Ne 


-3.91 


Ca 


-5.64 


Na 


-5.67 


Fe 


-4.33 


Mg 


-4.42 


Ni 


-5.75 



http:/ /www. sron.nl/spex 



To allow the cooling curve to be adapted for differ- 
ent abundances, we list the contributions of the differ- 
ent elements to the cooling rates, as calculated with the 
SPEX code. This is tabulated in Table [4] (online only). 
The cooling rates are gi ven again for solar abundances 
([Anders fc Grevessd 119891) . The contributions of the indi- 
vidual elements are plotted in Fig.[2]to illustrate the relative 
importance of the different elements to certain temperature 
ranges. The cooling rates are given for normalised values of 
nunc = 1, arid should be multiplied by n /riH when used 
in a hydrocode. 

For different metallicities, we can easily construct new 
cooling curves from Table[4] (online only), in which the cool- 
ing contributions and ionisation fractions of the elements 
are given separately. Cooling rates could then be calcu- 
lated by summing the contributions of the different ele- 
ments, multiplied by the appropriate abundance ratio and 
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Fig. 1. Cooling curves compared: The higher cooling rates 
calculated with SPEX are mainly due to a more complete 
coverage of the line transitions, including Fe L and EUV 
lines. 



ionisation fraction 

m 



rii (solar) 



A N (X h T), 



(3) 



where rij is the abundance ratio (relative to H) of ele- 
ment Xj. In this way, new cooling curves can be created 
ad hoc when needed for simulations of plasmas with differ- 
ent metallicities. 



3. Assumptions and validity of CIE cooling rates 

Collisional ionisation equilibrium is a valid assumption if 
an optically thin plasma is dominated by collisional pro- 
cesses and the cooling timescale is longer than either the 
ionisation or recombination timescales. Deviations from 
CIE can occur if there is additional photoionisation or 
there is non-equilibrium ionisation (NEI), i.e., the re- 
combination or ionisation timescales are longer than the 
cooling/heating timescale. These circumstances will pro- 
duce an overionised plasma if photoionisation is present 
or if the cooling timescale is shorter than the recombina- 
tion timescale. In young objects, such as supernova rem- 
nants, the heating timescale is shorter than the ionisation 
timescale, giving rise to an underionised plasma. 

NEI requires a time-dependent treatment of the cooling 
and is not easily implemented in standard hydrodynami- 
cal codes, but some specialised one-d imensional hydrody- 
namical codes for SN Rs do exist (e.g.. iBadenes et aLll2003t 
iKosenko et al.l [20041 ). On the other hand, photoionisation 
requires additional ionisation dependence on the photon- 
field, resulting in a higher ionisation than CIE for the same 
plasma temperature. The efficiency and reach of a photo- 
ionising source depends on both its spectrum and the ion- 
isation parameter £ = where R is the distance to the 
source. Depending on whether the elements are more or 
less effective coolants in their additionally ionised state, the 
cooling rates may be higher or lower, respectively. 

Studies of the cooling rates of overionised plasmas, rele- 
vant to photoionised plasmas or the late radiative stages 
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Fig. 2. Contributions of different elements to the cooling 
curve are given. Each of the plots shows a different set of 
elements. Important peaks are labelled with the name of 
the element. The total cooling curve (black solid line) is an 
addition of the individual elemental contributions. 



of supernova rem nants, have been performed by vari- 
ous authors (e.g.. iKafatod Il973t IShapiro fe Moord Il976t 
ISutherland fe Dopital 1993 ). The cooling rates in all of 
these cases are lower than that of a CIE cooling curve. 
For comparison, we plot the cooling curve by MB81 in 
Fig. [1] This curve u s es the time-dependent cooling curve 
by IShapiro fe Moord (| 19761 ) for the temperature range 4 < 
logT < 6. The gas is initially assumed to be in ionisa- 
tion equilibrium at T — 10 6 K, and cools time-dependently 
from there. Since the cooling timescale is shorter than the 
recombination timescale, the gas becomes overionised, lead- 
ing to suppressed cooling rates in this regime. Similarly, in 
intergalactic plasmas where photoio nisation is important , 
cooling rates tend to be suppressed (jWiersma et al.l l2009). 
Alternatively, cooling rates may be much higher when a 
neutral gas i s shock-heated, due to the initial ionisation 
of H and He (|Cox fe Raymond! 1 1985( 1. Unless full radiative 
transfer is taken into account, we therefore find CIE to be 
a reasonable assumption to use in circumstellar environ- 
ments and other low density plasmas. Only when photo- 
ionisation is significant, which depends on the tempera- 
ture of the star and the distance to the source, or when 
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strong shocks with velocities higher than about 150 km s^ 1 
are present, deviation from CIE can become important. 
At higher shock speeds, the gas is heated to temperatures 
higher than logT w 5.7, at which cooling by H and He 
lines of C, N, and O, becomes suppressed. Even in our rel- 
atively strong shock between two subsequent wind-phases 
as presented in Sect. [5] however, the typical speed at the 
forward shock is less than 150 km s" 1 and the temperature 
is typically around 10 5 K. 

Since recombination, collisional ionisation, and cooling 
have similar density dependences, the cooling rates that we 
calculate are independent of density, as long as the gas can 
be regarded as optically thin. 

4. Low temperature cooling 

The SPEX package is tailored to the X-ray regime, and 
includes continuum but no line emission for wavelengths 
longer than A > 2000 A. Therefore, it may underestimate 
cooling for temperatures below T = hc/Xk ~ 10 4 86 K. At 
temperatures below 10 4 K, cooling is much less efficient, 
and in the case of the typical CSM that we are interested 
in, most cooling takes place above this threshold. 

For low temperatures (T < 10 4 K), we use the 
cooling rates from a num ber of reactions as given by 
iDalgarno fc McCravl (|l972t ). In this regime, fractional ioni- 
sation and molecular line transitions can play an important 
role. We can set the ionisation fraction fi — n c /nn, which 
is a free parameter, to accomodate for different circum- 
stances in the circumstellar medium (CSM). This fraction 
strongly depends on specific conditions concerning shocks, 
radiation, and the history of the plasma. The cooling rates 
are included in Table [3] (online only) , and below we repeat 
the equations that we use to calculate them. The values for 
an ionisation fr action of 10~ 4 corr espond most closely to 
the values that ISmith et al.l ([2008) find for a Solar abun- 
dance plasma. 

The elements that are taken into account for low- 
temperature cooling are O, C, N, Si, Fe, Ne, and S. The 
cooling is provided partly by excitation of singly charged 
ions with the thermal electrons. The electron-collision 
transitions that are taken into account are given by the 
following equations (with the temperature T in Kelvin) 
(|Dalgarno fc McCravlfl972T) : 



A C (C+) 



A e (Si+) 



fAT- 1 ' 2 7.9 x 10- 20 e~ 92 / T 



n H 

3.0 X l -17 e -61900/T 



fi T- 1 ' 2 1.9 x 10- 18 fT 413 / T 



A e (Fe+) 



3.0 x iQ-^eSseoo/T 



^ fi T- 1 ' 2 4.8 x 10- ^e^/T 
n H 



1.1 x 10 -18 I e 



(4) 



(5) 



(6) 



r 554 / T + 1.3e- 961 / T ) 
+ 7.8 x i - 18 e 3496 / T 



Ae(0) 



fi) 1.74 x iq- 24 T 1/2 



[(1 - 7.6T- 1 / 2 )tr 228 / T 



(7) 



+ 0.38(1 - 7.7T- 1/2 )e~ 326/T ] 

+ 9.4 X l(r23 T l/2 e -22700/T 

A c (0+) = (j^f^ 1-5 x lQ-^T-^e" 38600 / 7 



A e (N) = — fi 8.2 x 10- 22 T 4 / 2 



x (1 - 2.7 x 10~ 9 T 2 )e 

n. 



-27700/T 



(8) 
(9) 

(10) 
(11) 

For low fractional ionisation, collisions with neutral hy- 
drogen can add substantially to the cooling. The following 
transitions are taken into account: 



Ac (S+) = ( ^fc ) 8.4 x io-ist- 1 /^- 21400 / 71 . 



A H (Si+) 
A H (Fe+) 



rij 

"H 

rij 
n H 

-554/T 



7.4 x i -23T- 1 / 2 e- 413 / T 
1.1 x 10~ 22 x 



(12) 
(13) 



1.4e 



-961/T 



The total cooling that we take into account for T < 10 4 K 
is a summation of contributions from Eq. 0] to 1131 plus 
contributions from collisions between electrons and hydro- 
en at oms as tabulated in Table 2 in IDalgarno fc McCravl 
1971) . 



5. Application: Test problem 

We use the AMRVAC code with the new cooling prescrip- 
tion to simulate the circumstellar medium (CSM) of a mas- 
sive star. The adaptive mesh refinement (AMR) ensures a 
computationally efficient approach to problems that have 
localised needs of high resolution, but where lower resolu- 
tion suffices in substantial parts of the grid. The standard 
Euler equations are solv ed in a con servative way, where we 
apply an HLLC solver (|Torol [T999h to capture more accu- 
rately the contact discontinuity. 

We implement the stages of the evolution of a mas- 
sive star in which a red supergiant (RSG) wind is fol- 
lowed by the much faster Wolf-Rayet (WR) wind. A RSG 
wind is typically slow, cool, and dense because of the 
low escape velocity from the surface of the giant star. 
The WR wind, however, is fast and hot, originating from 
the compact stellar core. The WR wind therefore sweeps 
up the RSG wind into a shell (|Garcia-Segura fc Mac Low! 
119951 : iGarcia-Segura et all 119961: iDwarkadasl 12007ft . here- 
after called the WR shell, whose thickness and stability 
is strongly affected by the inclusion of radiative cooling. 

We set up the simulation in 2.5D on a spherical grid 
in a meridional plane over a full angle 7T, using reflective 
boundary conditions and symmetry around the 0-axis. The 
radial extent is about 3 pc and we use a base resolution 
of 60 x 48 gridcells with 4-6 refinement levels. For each 
refinement level, the resolution is doubled. The values that 
we use for the wind parameters are the following: Mrsg — 



"\ T RSG - 10 
1.7 x 10 3 km s 



3 K, 
-l 



1.54 x 10~ 5 Af o yr _1 , v RS g = 4.7 km s 
M W r = 9.7 x 1O- 6 M yr" 1 , v WR = 
T RSG = 10 4 K. 

For the RSG phase, 4 refinement levels are used, and 
for the WR phase, we use 4-7 refinement levels and check 
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for differences in the development of the instabilities and 
the resolution dependence of the radiative losses. The max- 
imum resolution corresponds to 1.56 x 10 15 cm x 0.0156°. 

In general, all plasmas are prone to radiative losses. 
Depending on the circumstances however, these can induce 
significant changes in the hydrodynamics, or have only mi- 
nor quantitative effects. In the situation that we simulate, 
without cooling the maximum compression reached in the 
shell is a factor of four; the value expected for a strong shock 
in a mono-atomic gas with an adiabatic index of 7 = 5/3. 
With cooling, the shell can become much more compressed 
and, as a result, much thinner. When the compression ra- 
tio exceeds a factor of about 21, the shell c an become sub- 
ject to the Vishniac thin shell instability (|Vishnia3 Il983t 
iRvu fc Vishniadll987trMac Low k Normanll 19931 ). 

The Vishniac thin shell instability, which actually is an 
overstability, occurs when a thin shell is perturbed while 
being driven by a high pressure region from the inside and 
encountering ram pressure on the outside. The high pres- 
sure from the WR wind acts in a perpendicular direction 
to the shell, whereas the ram pressure acts in a parallel di- 
rection to the velocity. Whereas the velocity is radial, the 
interior pressure in the case of a perturbed shell can create a 
flow along the shell surface, driving matter from advanced 
regions to regions lagging behind. The more tenuous re- 
gions then become more susceptible to deceleration by ram 
pressure, reversing the cycle. 

Above a temperature of 10 5 K, the plasma also becomes 
susceptible to the thermal instability. This occurs for shock 
veloc ities above 150 km s" 1 , when (din A)/ (din T) < 2 
(e^ lFieldl Il965t lBertschingerlll986l : ISutherland fc Dopital 
1993). For typical circumstellar wind parameters, the 
plasma temperature in the cooling shell remains below this 
regime. 

The cooling of the plasma is calculated at each timestep 
and the thermal energy is correspondingly updated. The 
temperature of the plasma is first determined. The corre- 
sponding cooling rate is then found by interpolating be- 
tween the appropriate values in the cooling curve. The en- 
ergy loss rate per unit volume is calculated as follows: 



de 
dt 



= -A hd nH- 



(14) 



The energy losses are subtracted from the thermal energy 
only, preventing numerical precision problems from produc- 
ing negative energy. 

Although the cooling rates are higher when we use the 
SPEX-calculated cooling curve, in both cases the thickness 
of the WR shell is not so much limited by the cooling curve 
used, as by the adopted numerical resolution. It is impor- 
tant however, that the shell becomes affected by instabili- 
ties in the presence of cooling, whereas it remains stable in 
the absence of cooling. We compare the differences for runs 
with different resolution and different cooling prescription. 
For comparison, we use the cooling curve as prescribed in 
IMacDonald k Bailevl (|1981l ). since this is the cooling curve 
that is used in the ZEUS3D v3.4 code, as used by some 
of the authors befo re, for simulations of a WR wind in 
iSchure et al.1 (f2008h . 

Figure [3] shows the density of the circumstellar medium 
after the WR-phase has just set in for about 1750 years, and 
shows the early developments of the Vishniac thin shell in- 
stability. The fast wind from the Wolf-Rayet star has swept 
up the red supergiant wind of the previous evolutionary 




logrtio 



n.4 n.fi 11. 



Fig. 3. WR wind colliding with RSG wind. The density 
of the CSM is shown about 1750 years after the initia- 
tion of the WR wind, showing the early development of 
the Vishniac thin shell instability. The left panel shows the 
CSM when no radiative losses are taken into account. The 
middle panel shows the case where the MB81 cooling pre- 
scription is applied, and the right panel shows the CSM 
with the cooling prescription as presented in this work. The 
black line indicates the contact discontinuity. The num- 
ber of refinement levels used for the resolution is 5, i. e., 
6.25 x 10 15 cm x 0.0625°. 



stage and four regions can be identified. From inside (R=0) 
out: the free streaming fast WR wind, the shocked WR 
wind. The black line indicates the contact discontinuity, 
behind which the shocked, and farther out the free stream- 
ing, RSG wind can be found. Only part of the grid is shown 
to enlarge the interesting region. 

The left panel shows the CSM for the colliding winds 
in the absence of cooling. The shell is smooth and thick, 
with a maximum compression ratio that is similar to that 
of a strong shock, i. e., a factor 4. The forward shock is 
partly outside the plotted region of the grid and extends to 
a radius of about 1.2 x 10 18 cm. The middle and right pan- 
els illustrate r esults when cooling is present , for the cooling 
rate given by IMacDonald fc Bailevl (|l98ll ) in the middle 
panel, and for the cooling rate calculated with the SPEX 
code in the right panel. The instability is induced by ran- 
dom density perturbations of the 1% level, which may in- 
duce little differences between runs. If we do not include the 
initial perturbations, the code retains the original spheri- 
cal symmetry to high accuracy, and the instabilities do not 
arise until very late. The extrapolation of the MB81 cooling 
curve towards lower temperatures results in higher cooling 
rates below 10 4 K. This gives rise to more cooling in the 
compressed RSG shell and a more patchy structure as can 
be seen in the middle panel of Fig. [3] For temperatures 
above 10 4 K, cooling rates are higher for the SPEX curve. 
This applies to the region of shocked WR wind material 
and results in this area being more compressed in the right 
panel. Consequently, the WR shell in the middle panel ex- 
tends a little farther out. 

Total energy losses are higher with our new cooling pre- 
scription. The comparison is made in Fig. 2] The thermal 
energy per unit volume is plotted for the different cases in 
Fig. [5] The inclusion of cooling reduces the thermal energy 
substantially. We note again that the initial random pertur- 
bation will influence the later evolution, explaining some of 
the inter-resolution differences. However, the key point is 
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the systematic reduction between the SPEX versus MB81 
cooling prescriptions. 



no cooling 
MB81 cooling < 
this work: nl 5 
this work: nl 4 
this work: nl 6 
this work: nl 7 




0.5 - 



0.0 



1000 
time (yr) 



1500 



Fig. 4. Thermal energy in the stellar winds. The black solid 
line shows the energy when no cooling is taken into account. 
The black dotted curve shows the energy when the MB81 
cooling rates are used, and the black dashed curve shows 
the system for the cooling rates presented in this work. The 
gray lines show the resolution dependence of the thermal 
energy, where nl indicates the number of refinement levels 
used. 

The wavelength of the Vishniac thin shell instability de- 
pends on the thickness of the shell. In Fig. we show the 
effect of resolution on the development of the instabilities 
with the new cooling curve. For higher resolution, the shell 
becomes thinner, and as a result the wavelength of the in- 
stability shortens. Although the WR shell is not resolved 
in terms of the dynamics, we see that for thermal losses, 
convergence seems to be reached for the resolutions used. 
Figure [4] shows the resolution dependence of the thermal 
losses for the different cooling prescriptions and different 
resolutions. 

A typical WR phase lasts for about 10 5 years. In Fig. [6j 
we show the development of the WR shell later in the evo- 
lutionary stage, at which time the WR phase has lasted 
for about 22,700 yr. It is clear that the presence of cool- 
ing majorly alters the CSM evolution, by losing energy by 
radiation, and by the irregular shell that develops as a re- 
sult of this. At later times, shells from earlier evolutionary 
phases may intervene with the CSM evolution and simu- 
lations need to t ake these phases into account (see e.g., 
lDwarkadasll2007Hvan Marie et al.ll2007l ). 

6. Conclusion and discussion 

SPEX is a very complete package that calculates the emis- 
sivity of a plasma at a certain temperature, taking into 
account the abundances. We calculate the radiative losses 
for a plasma in CIE, and com pare these with the curve 
bv lSutherland fc Dopital (|l993l ). which was calculated with 
the same basic assumptions. We find that cooling rates 
calculated by SPEX tend to be higher than those of 
ISutherland fc Dopital (|l993t ). This is not surprising, since 
more lines have since been added and the SPEX package is, 




0.2 0.4 0.6 0.H 1 0.2 0.4 0.6 O.fi 1 



0.2 IL4 11.6 II. X 



Fig. 5. The Vishniac thin shell instability depends on nu- 
merical resolution. At high resolution, the shell becomes 
thinner and thus the wavelength shorter. The different pan- 
els correspond to refinement levels of 4, 6, and 7, corre- 
sponding to maximum resolution of 1.25 x 10 16 cm x 0.125°, 
3.125 x 10 15 cm x 0.03125°, and 1.56 x 10 15 cm x 0.0156°. 
The simulation with 5 refinement levels is shown in Fig. [31 



logrho 




2 4 6 8 
R (10*18 cm) 



Fig. 6. Unstable WR shell at a time of about 22,700 yrs 
after its transition to a WR star. The shell becomes highly 
irregular due to the Vishniac thin shell instability. The 
black line indicates the contact discontinuity. This simu- 
lation used 6 refinement levels for the resolution. 



to our knowledge, the most complete package for X-rays and 
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EUV. F or low temperatures, we reu se the cooling prescrip- 
tion by iDalgarno fc McCravl (|1972f ). Using that, and the 
SPEX results, we provide a recommended cooling curve for 
Solar abundances^!, which can be adjusted to different ioni- 
sation fractions at low temperatures. Additionally, changes 
can be decided for each cooling curve according to spe- 
cific requirements on the abundances, by using the separate 
cooling rates of the individual elements. 

In our application, we calculate the CSM of a massive 
star that experiences a RSG phase followed by a WR phase. 
The colliding winds form a shell that, in the presence of 
cooling, becomes very thin and susceptible to the Vishniac 
thin shell instability. Taking into account the cooling of 
the plasma is very important when accurately simulating 
systems such as these. 
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Table 2. Cooling curve for solar metallicity 



log T log A N log A hd n e /n H 

(K) (erg s _1 cm 3 ) (erg s _1 cm 3 ) 



Wiersma, R. P. C, Schaye, J., & Smith, B. D. 2009, MNRAS, 393, 99 



3.80 


-25.7331 


-30.6104 


1.3264E-05 


3.84 


-25.0383 


-29.4107 


4.2428E-05 


3.88 


-24.4059 


-28.4601 


8.8276E-05 


3.92 


-23.8288 


-27.5743 


1.7967E-04 


3.96 


-23.3027 


-26.3766 


8.4362E-04 


4.00 


-22.8242 


-25.2890 


3.4295E-03 


4.04 


-22.3917 


-24.2684 


1.3283E-02 


4.08 


-22.0067 


-23.3834 


4.2008E-02 


4.12 


-21.6818 


-22.5977 


1.2138E-01 


4.16 


-21.4529 


-21.9689 


3.0481E-01 


4.20 


-21.3246 


-21.5972 


5.3386E-01 


4.24 


-21.3459 


-21.4615 


7.6622E-01 


4.28 


-21.4305 


-21.4789 


8.9459E-01 


4.32 


-21.5293 


-21.5497 


9.5414E-01 


4.36 


-21.6138 


-21.6211 


9.8342E-01 


4.40 


-21.6615 


-21.6595 


1.0046 


4.44 


-21.6551 


-21.6426 


1.0291 


4.48 


-21.5919 


-21.5688 


1.0547 


4.52 


-21.5092 


-21.4771 


1.0767 


4.56 


-21.4124 


-21.3755 


1.0888 


4.60 


-21.3085 


-21.2693 


1.0945 


4.64 


-21.2047 


-21.1644 


1.0972 


4.68 


-21.1067 


-21.0658 


1.0988 


4.72 


-21.0194 


-20.9778 


1.1004 


4.76 


-20.9413 


-20.8986 


1.1034 


4.80 


-20.8735 


-20.8281 


1.1102 


4.84 


-20.8205 


-20.7700 


1.1233 


4.88 


-20.7805 


-20.7223 


1.1433 


4.92 


-20.7547 


-20.6888 


1.1638 


4.96 


-20.7455 


-20.6739 


1.1791 


5.00 


-20.7565 


-20.6815 


1.1885 
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-20.7820 
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5.68 


-21.5098 
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-21.7696 


-21.6886 


1.2050 



2 http:/ /www. phys.uu.nl/~schure/cooling 
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Table 3. Cooling curve for solar metallicity fo r T < 10 4 K, 
as computed from Dalgarno & McC ravl (|1972t ). see Sect. 2] 
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Table 4. Cooling rates per clement for solar metallicity 
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,04 


-23. 


,46 


-24.29 


-24. 


.45 


-22 


,02 


-24.66 


5 


.96 


-23. 


.58 


-23 


.56 


-23 


.56 


-24 


.32 


-23 


.33 


-23 


.07 


-24, 


,39 


-22.48 


-23 


,94 


-22. 


,96 


-23. 


,44 


-24.22 


-24. 


,52 


-22. 


,01 


-24.48 



o 



T H He C N O Ne Na Mg Al Si S Ar Ca Fc Ni 

6.00 -23.57 -23.58 -23.52 -24.23 -23.37 -23.26 -24.61 -22.57 -24.00 -22.88 -23.42 -24.17 -24.51 -22.01 -24.32 

6.04 -23.57 -23.60 -23.49 -24.14 -23.36 -23.42 -24.82 -22.74 -24.05 -22.82 -23.38 -24.15 -24.46 -22.01 -24.17 

6.08 -23.56 -23.62 -23.50 -24.05 -23.33 -23.56 -25.01 -22.95 -24.11 -22.77 -23.35 -24.14 -24.40 -22.02 -24.04 

6.12 -23.55 -23.64 -23.53 -23.98 -23.26 -23.68 -25.18 -23.17 -24.19 -22.74 -23.32 -24.15 -24.33 -22.02 -23.92 

6.16 -23.54 -23.65 -23.58 -23.93 -23.19 -23.79 -25.33 -23.37 -24.31 -22.73 -23.30 -24.17 -24.28 -22.03 -23.82 

6.20 -23.53 -23.66 -23.64 -23.90 -23.11 -23.86 -25.47 -23.56 -24.47 -22.77 -23.29 -24.20 -24.25 -22.05 -23.74 

6.24 -23.52 -23.67 -23.71 -23.89 -23.04 -23.91 -25.59 -23.72 -24.64 -22.87 -23.29 -24.24 -24.24 -22.09 -23.67 

6.28 -23.51 -23.68 -23.79 -23.91 -22.98 -23.94 -25.68 -23.86 -24.81 -23.02 -23.29 -24.28 -24.23 -22.17 -23.63 

6.32 -23.50 -23.68 -23.86 -23.95 -22.93 -23.93 -25.75 -23.98 -24.96 -23.19 -23.29 -24.31 -24.23 -22.29 -23.61 

6.36 -23.49 -23.68 -23.93 -24.01 -22.91 -23.90 -25.78 -24.08 -25.11 -23.36 -23.32 -24.32 -24.22 -22.44 -23.63 

6.40 -23.47 -23.69 -24.00 -24.07 -22.91 -23.86 -25.78 -24.16 -25.23 -23.52 -23.39 -24.32 -24.19 -22.58 -23.72 

6.44 -23.46 -23.69 -24.06 -24.14 -22.94 -23.80 -25.76 -24.22 -25.33 -23.67 -23.50 -24.31 -24.16 -22.69 -23.87 

6.48 -23.45 -23.68 -24.12 -24.21 -22.98 -23.75 -25.71 -24.25 -25.41 -23.79 -23.63 -24.31 -24.14 -22.76 -23.98 

6.52 -23.43 -23.68 -24.18 -24.28 -23.04 -23.70 -25.65 -24.26 -25.47 -23.90 -23.76 -24.35 -24.12 -22.80 -24.04 

6.56 -23.42 -23.68 -24.23 -24.34 -23.10 -23.66 -25.58 -24.24 -25.49 -23.98 -23.89 -24.43 -24.12 -22.81 -24.05 

6.60 -23.40 -23.67 -24.28 -24.41 -23.17 -23.64 -25.52 -24.22 -25.49 -24.05 -24.00 -24.55 -24.14 -22.80 -24.05 

6.64 -23.39 -23.67 -24.32 -24.46 -23.23 -23.63 -25.47 -24.18 -25.46 -24.09 -24.10 -24.68 -24.20 -22.78 -24.04 

6.68 -23.37 -23.66 -24.36 -24.52 -23.30 -23.64 -25.42 -24.14 -25.42 -24.12 -24.19 -24.81 -24.29 -22.74 -24.02 

6.72 -23.36 -23.65 -24.40 -24.57 -23.36 -23.66 -25.39 -24.10 -25.37 -24.13 -24.26 -24.91 -24.41 -22.71 -24.02 

6.76 -23.34 -23.64 -24.44 -24.61 -23.41 -23.70 -25.38 -24.07 -25.31 -24.12 -24.31 -25.00 -24.55 -22.68 -24.03 

6.80 -23.33 -23.64 -24.47 -24.66 -23.47 -23.75 -25.39 -24.05 -25.26 -24.10 -24.34 -25.06 -24.67 -22.65 -24.06 

6.84 -23.31 -23.63 -24.50 -24.70 -23.52 -23.80 -25.41 -24.03 -25.21 -24.07 -24.36 -25.11 -24.79 -22.64 -24.10 

6.88 -23.29 -23.61 -24.53 -24.73 -23.56 -23.86 -25.44 -24.04 -25.17 -24.05 -24.37 -25.14 -24.89 -22.64 -24.15 

6.92 -23.28 -23.60 -24.55 -24.77 -23.61 -23.92 -25.48 -24.05 -25.15 -24.02 -24.36 -25.15 -24.97 -22.65 -24.20 

6.96 -23.26 -23.59 -24.57 -24.80 -23.65 -23.97 -25.53 -24.08 -25.14 -24.00 -24.35 -25.15 -25.04 -22.67 -24.25 

7.00 -23.24 -23.58 -24.59 -24.82 -23.68 -24.02 -25.58 -24.12 -25.14 -23.98 -24.33 -25.13 -25.09 -22.70 -24.29 

7.04 -23.23 -23.57 -24.61 -24.85 -23.72 -24.07 -25.63 -24.16 -25.15 -23.98 -24.30 -25.11 -25.12 -22.74 -24.32 

7.08 -23.21 -23.55 -24.62 -24.87 -23.75 -24.12 -25.68 -24.21 -25.18 -23.98 -24.28 -25.07 -25.14 -22.79 -24.35 

7.12 -23.19 -23.54 -24.63 -24.89 -23.77 -24.16 -25.73 -24.26 -25.22 -24.00 -24.26 -25.04 -25.15 -22.86 -24.39 

7.16 -23.17 -23.52 -24.64 -24.91 -23.80 -24.20 -25.77 -24.30 -25.26 -24.02 -24.25 -25.00 -25.15 -22.94 -24.45 

7.20 -23.16 -23.51 -24.65 -24.92 -23.82 -24.23 -25.81 -24.35 -25.30 -24.06 -24.24 -24.96 -25.13 -23.03 -24.51 

7.24 -23.14 -23.50 -24.66 -24.94 -23.84 -24.27 -25.85 -24.40 -25.34 -24.09 -24.24 -24.93 -25.11 -23.13 -24.58 

7.28 -23.12 -23.48 -24.66 -24.95 -23.86 -24.30 -25.89 -24.44 -25.39 -24.13 -24.25 -24.90 -25.09 -23.21 -24.64 

7.32 -23.10 -23.46 -24.67 -24.96 -23.87 -24.32 -25.92 -24.47 -25.43 -24.18 -24.27 -24.88 -25.07 -23.29 -24.70 

7.36 -23.08 -23.45 -24.67 -24.97 -23.88 -24.35 -25.95 -24.51 -25.47 -24.22 -24.30 -24.87 -25.05 -23.36 -24.75 

7.40 -23.07 -23.43 -24.67 -24.97 -23.90 -24.37 -25.97 -24.54 -25.50 -24.26 -24.33 -24.86 -25.03 -23.42 -24.80 

7.44 -23.05 -23.42 -24.67 -24.98 -23.90 -24.39 -26.00 -24.57 -25.54 -24.29 -24.36 -24.87 -25.01 -23.46 -24.84 

7.48 -23.03 -23.40 -24.67 -24.98 -23.91 -24.41 -26.02 -24.60 -25.57 -24.33 -24.40 -24.88 -25.00 -23.49 -24.87 

7.52 -23.01 -23.38 -24.66 -24.98 -23.92 -24.42 -26.04 -24.62 -25.60 -24.36 -24.43 -24.90 -25.00 -23.52 -24.90 

7.56 -22.99 -23.37 -24.66 -24.98 -23.92 -24.44 -26.05 -24.64 -25.62 -24.39 -24.47 -24.92 -25.00 -23.53 -24.92 

7.60 -22.97 -23.35 -24.65 -24.98 -23.92 -24.45 -26.07 -24.66 -25.65 -24.42 -24.50 -24.95 -25.01 -23.54 -24.93 

7.64 -22.95 -23.33 -24.65 -24.97 -23.93 -24.45 -26.08 -24.68 -25.67 -24.45 -24.53 -24.98 -25.02 -23.54 -24.94 

7.68 -22.94 -23.32 -24.64 -24.97 -23.92 -24.46 -26.09 -24.70 -25.68 -24.47 -24.56 -25.01 -25.04 -23.54 -24.94 

7.72 -22.92 -23.30 -24.63 -24.97 -23.92 -24.47 -26.10 -24.71 -25.70 -24.49 -24.59 -25.03 -25.06 -23.54 -24.94 

7.76 -22.90 -23.28 -24.62 -24.96 -23.92 -24.47 -26.10 -24.72 -25.71 -24.51 -24.61 -25.06 -25.08 -23.53 -24.93 

7.80 -22.88 -23.26 -24.61 -24.95 -23.92 -24.47 -26.11 -24.73 -25.72 -24.52 -24.64 -25.09 -25.10 -23.53 -24.93 

7.84 -22.86 -23.25 -24.60 -24.95 -23.91 -24.47 -26.11 -24.73 -25.73 -24.53 -24.66 -25.11 -25.12 -23.53 -24.92 

7.88 -22.84 -23.23 -24.59 -24.94 -23.91 -24.47 -26.11 -24.74 -25.74 -24.55 -24.67 -25.13 -25.15 -23.52 -24.91 

7.92 -22.82 -23.21 -24.58 -24.93 -23.90 -24.47 -26.12 -24.74 -25.75 -24.56 -24.69 -25.15 -25.17 -23.52 -24.90 

7.96 -22.80 -23.19 -24.57 -24.92 -23.89 -24.47 -26.11 -24.75 -25.75 -24.56 -24.70 -25.17 -25.19 -23.53 -24.89 

8.00 -22.78 -23.17 -24.56 -24.91 -23.88 -24.46 -26.11 -24.75 -25.76 -24.57 -24.71 -25.19 -25.21 -23.53 -24.89 

8.04 -22.76 -23.16 -24.54 -24.90 -23.87 -24.46 -26.11 -24.75 -25.76 -24.57 -24.72 -25.20 -25.22 -23.54 -24.89 

8.08 -22.74 -23.14 -24.53 -24.88 -23.86 -24.45 -26.10 -24.74 -25.76 -24.57 -24.73 -25.21 -25.24 -23.55 -24.89 

8.12 -22.72 -23.12 -24.52 -24.87 -23.85 -24.45 -26.10 -24.74 -25.75 -24.58 -24.73 -25.22 -25.25 -23.56 -24.89 

8.16 -22.70 -23.10 -24.50 -24.86 -23.84 -24.44 -26.09 -24.74 -25.75 -24.58 -24.74 -25.23 -25.27 -23.57 -24.89 



